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We study the scattering problem for the non-Hermitian scattering center, which consists of two 
Hermitian clusters with anti-Hermitian couplings between them. Counterintuitively, it is shown that 
it acts as a Hermitian scattering center, satisfying \r\ 2 + \t\ 2 = 1, i.e., the Dirac probability current is 
conserved, when one of two clusters is embedded in the waveguides. This conclusion can be applied 
to an arbitrary parity-symmetric real Hermitian graph with additional PT-symmetric potentials, 
which is more feasible in experiment. Exactly solvable model is presented to illustrate the theory. 
Bethe ansatz solution indicates that the transmission spectrum of such a cluster displays peculiar 
feature arising from the non-Hermiticity of the scattering center. 
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I. INTRODUCTION 

A non-Hcrmitian Hamiltonian is usually endowed with 
the physical meaning when it possesses entirely real quan- 
tum mechanical energy spectrum and the complex exten- 
sion of the conventional quantum mechanics, a parity- 
time (TT) symmetric quantum theory, has been well de- 
veloped [l|-|8( since the seminal discovery by Bender [if. 
Such a theory gives the pseudo-Hermitian Hamiltonian a 
physical meaning via its corresponding Hermitian coun- 
terparts [9Hllj|. which has an identical spectrum. The 
metric-operator theory outlined in Ref. 0] provides a 
mapping of such a pseudo-Hermitian Hamiltonian to an 
equivalent Hermitian Hamiltonian. Thus, most of the 
studies focused on the quasi-Hermitian system, or un- 
broken 'PT-symmetric region [l2|, [l3j]. However, the ob- 
tained equivalent Hermitian Hamiltonian is usually quite 
complicated [(| [l4| , involving long-range or nonlocal in- 
teractions, which is hardly realized in practice. 

Experimentally, the PT symmetry is of great relevance 
to the technological applications based on the fact that 
the imaginary potential could be realized by complex in- 
dex in optics |15l - [l8j . Furthermore, the PT optical po- 
tentials can be realized through a judicious inclusion of 
index guiding and gain/loss regions. Such non-Hermitian 
systems are not isolated but usually embedded in the 
large Hermitian waveguides. Pure imaginary potential 
as a scattering center breaks the conservation of the flow 
of probability Thus, it is interesting to investigate 
what happens when the non-Hcrmitian system is with 
balanced gain and loss as a scattering center, and much 
effort devoted to such a topic is based on the framework 
of rj- metric [2CH23| . 

In this paper, we study the scattering problem for the 
non-Hcrmitian scattering center based on the configura- 
tions involving two arbitrary Hermitian networks cou- 
pled with anti-Hermitian interaction. It is shown that 
for any scattering state of such a non-Hermitian system, 
the Dirac probability current is always conserved at any 



degree of the non-Hermiticity. We apply such a rigor- 
ous result to the system with PT-symmetric potentials, 
which is more feasible in experiment. 

This paper is organized as follows. Section [TTI presents 
the exact analytical solution of the scattering problem for 
the concerned non-Hermitian scattering center. Section 
IIIII is the application of the rigorous result to the system 
with PT-symmetric potentials. Section IIVI consists of 
an exactly solvable example to illustrate our main idea. 
Section [V] is the summary and discussion. 



II. MODEL AND FORMALISM 

In general, a non-Hermitian Hamiltonian H is related 
by a similarity transformation to an equivalent Hermi- 
tian Hamiltonian h. Such a connection is valid within 
the so called unbroken symmetric region. However, when 
a non-Hermitian system interacts with other Hermitian 
system, such a region loses its physical meaning: On the 
one hand, the unbroken symmetric region is shifted in the 
whole non-Hermitian system. On the other hand, it may 
act as a Hermitian system in the scattering problem with- 
out the restriction on the degree of the non-Hermiticity. 
In this section, we will investigate the latter situation. 

The Hamiltonian of the concerned scattering tight- 
binding network has the form 



H — Hl + Hr + He , 



(1) 



where 



H l = - k Y, \3)l 0* — 11 — 9L \L) L + H.c, (2) 

+ 00 

Hr = \j) R + 9R \R) R (1| + H.c, (3) 

represent the left (Hl) and right (Hr) waveguides with 
real k and 
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He = Ha + Hb + Hab + Hba , 



(4) 
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FIG. 1. ((Color online) Schematic illustration of the config- 
uration of the concerned network. It consists of two arbi- 
trary graphs of the Hermitian tight-binding networks A and 
B (shadow) with one of them connecting to two semi-infinite 
chains as the waveguides at the joint sites L and R. The non- 
Hermiticity of the whole scattering center A + B arises from 
the anti-Hermitian interaction (dished lines) between them. 
It is shown that the non-Hermitian scattering center acts as 
a Hermitian one, preserving the Dirac probability current. 



describes a non-Hermitian network as a scattering center. 
Here \L) and \R) denote the sites state on the joint sites 
on the network A, which are simply taken as \L) = |1) * 
and \R) — \Na) a without losing the generality. The 
subgraphs 



N A 

H A = (HA) t3 \i) A 01 , 

N B 

h b = 2 (H B ) tJ \i) B oi , 



(5) 
(6) 



are arbitrary Hermitian networks, i.e., H A = H A , and 
H B = Hb, while the coupling between them is anti- 
Hermitian, i.e., H AB = —Hba- 

Na N b 

H AB=J2Y,( H * B k\ i )AB(j\i (7) 
i=l j=l 

then the scattering center with respect to the basis 
{\i) A , \i) b\ IS m ^ ne form of 



— oo Na N b +oo 

M - £ U li>z+I> Ij>a+X> \i) B +Efj li> 

i=-i i=i i=i j=i 



R ■ 



where the scattering wavefunction fj is in form of 



+ re 



-ikj 



te ifci ,0' > 1) 



(9) 



(10) 



Here r, i are the reflection and transmission coefficients 
of the incident wave, which is what we concern only in 
this paper. Substituting the wavefunction \ijjk) into the 
Schrodinger equation 



H\1> k ) = E\j, k ), 



(11) 



the explicit form of the Schrodinger equations in 
the truncated Hilbert space spanned by the basis 
G [l,N A }) A ,\j,j e [1,N B }) B } can be expressed in 
the following matrix equation form 



^Na 
Pi 



/ 9hf-i \ 









V o J 



(12) 



where A is an (Na + Nb) X (JVa + 2Vb) matrix defined 
by 



A 



Ha — E H AB 
~ h ab h b- E 



(13) 



From the reduced Schrodinger equation of Eq. (|12p , we 
obtain 



, H A H ab 

H C = I rrt 



(8) 



The non-Hermiticity of He arises from this anti- 
Hermitian term. The non-Hermitian Hamiltonian He 
may have fully real spectrum or not. In the following, we 
will show that it always acts as a Hermitian scattering 
center no matter the reality of the spectrum. 

For an incident plane wave with momentum k incoming 
from the left waveguide L with energy E = —2k cos k, 
the scattering wave function \tpk) can be obtained by the 
Bethe ansatz method. The wave function has the form 



a x = (A- 1 ) n g L f_ 1 + (A' 1 ) g R fi, (14) 



a N . 



n a n a 



9B.fl, 



Here A 1 is the inverse of matrix A, with the element 
being expressed as 



(A- 1 ), 



Cji 



det (A) 



(-1)' +J det (Mjj) 
det (A) 



(15) 



in term of the matrix of cofactors Cy. Here is the 
matrix obtained by deleting the ith row and jth column 
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from the matrix A. On the other hand, the Schrodinger 
equations for the sites of the waveguides connecting to 
the joints of the scattering center are 



- «/-2 - 9*L a l = Ef-x, 
-K.f 2 - 9*rO-n a = Efx, 



(16) 



which lead to 



Q i = — C 1 + r ) > a ^A 
9l 



9r 



Then associating with Eqs. (fl4|) . we have 



r = (-bb + ac — ae lk - ce l 
t = i2bsink/r], 



(17) 



(18) 
(19) 



where 



?7 = (66 - ac)e l2k + {a + c) e ik - 1, 

a = (A-% \g L \ 2 /«, c = (&- 1 ) NaNa \ 9 r\ 2 /«, (20) 

*=(*~ 1 )iN A Sl9R/K,b= {^ 1 ) NAl g L9 * R /n. 

One can determine the unknown coefhcients a, b, b, c 
and 77 through the matrix A by requiring that invertible 
matrix (Ha — E) or (Hb — E) exists. 
In the Appendix, we will show that 



fA- 1 ) 



A" 



for i,j £ [1, Na], or more explicitly for special cases 



(21) 



fA- 1 ) 



A" 



' 11 ' 



(A" 1 



N A N A I" >N A N A ' 



IN A 



(^~)n a i ' 



(22a) 
(22b) 
(22c) 



which indicate that both a and c are real, and b — b*. It 
is somewhat surprising that we get the conclusion from 
Eqs. (p j) . (|T9 ]) , (|20 ]l and p2 j) that 



(23) 



1*1 =1, 



which is common phenomenon in a Hermitian system but 
surprising in a non-Hermitian system. 



III. PT-SYMMETRIC POTENTIALS 

The accessible setup of non-Hermitian system in the 
lab is the PT-symmetric potentials, which can be real- 
ized through a judicious inclusion of index guiding and 
gain/loss regions. In the following, we will apply the ob- 
tained result to the system with the PT-symmetric po- 
tentials, in which the PT-symmetrical axis is along the 
waveguides. 



The geometry of the scattering center contains N\ 
2N 2 sites and possesses the following symmetry, 



P:|j> c — >]J>c> 0'e[l,Wi]) 

V : |i) c — ► \j) c = \N 2 + j) c , (j - N! £ [1,N 2 ]) 



(24) 



with the joint points L, R £ [l,Ni], where \j) c is the 
mirror symmetric counterpart of state \j) . We define 
the Hamiltonian of the center has the form 



^1+2^2 

Hvt= E ^K> c 01+H.c. (25) 

»,i=i)(*<i) 

JVi JVi+JV 2 

+E c/ ,-ii) c oi + E ^\j) c (j\+y;\j)M), 

3 = 1 j=Nx+X 

where Ky and Uj are real. In the Hilbert space spanned 
by basis {\j) c } (j £ [1,N 1 +2N 2 ]), the matrix of the 
Hamiltonian H-pr has the form 



H~ a 



H vr = H\ a H a + H S H a 



H 



(26) 



H, 



a/3 



H a — Hs 



where 



(Hsh = SisW - V*)/2 = ilm (V,) 8 K 



(27) 



and H 1 (H a ) is an Ni (N 2 ) dimension square matrix. 
The matrices H y , H a , and H a p are all real Hermitian 
while H ia is real. We can see that Hamiltonian H-pf 
describes an arbitrary real Hermitian graph with parity- 
symmetry as defined in Eq. (|24p combining with the on- 
site 'PT-symmetric potentials Hs- Thus Hpj- satisfies 
[PT,#pt]=0. 

Introducing the linear transformation 



\J)a = 



(ieli,^]) 



(li>c + li> c )/V2, [j-Ni£[l,N 2 ]) 



(28a) 



\J) B = (li) c - \3) e )/V2, (j -N,£ [1,N 2 ]) (28b) 



one can rewrite the matrix of Eq. (1261) in the basis 
J € [l,Ni + N 2 ]) A , |i, j £ [m + l,m + N 2 ]) B } as 
the form 



/ H 7 V2H ia 
H VT = \y/2H\ a H a +H a(s H s \ 
\ Hs H a — H a p 

Obviously, it is the special case of Eq. ©, where 



_ ( H 1 \p2H la 
A ~\ V2H\ a H a + H aP 



Hb = H, 
Hab = 



H 



a/3, 







,Hba={0 H s ) 



(30) 
(31) 
(32) 
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FIG. 2. (Color online) Schematic illustration for the exem- 
plified system, (a) 4-site non-Hermitian scattering center con- 
figuration, which consists of two on-site imaginary potentials 
171 and —172 connecting to two semi-infinite chains as the 
waveguides at the joint sites |1} C and |3) c , with the hopping 
strength —k. (b) The equivalent Hamiltonian of He [Eq. 
(j48}], which is obtained under the linear transformation of 
Eq. (|47|) . The hopping strengths between the site \A) and 
|1) , 1 3) are both — \/2/t. The dashed (green) line represents 
the effective hopping between sites \A) and \B) which is pure 
imaginary i (71 + 72) /2, with both potentials on \ A) and \B) 
being i (71 — 72) /2. In the case of 71 =72, it is shown that 
the non-Hermitian scattering center acts as a Hermitian one, 
preserving the Dirac probability current. 



written as 

—00 

H l = -kJ2 U - 1| - « |-1)lc (11 + H.c, (37) 
j=-i 

H R = -«X>*>* (j + 1| - k \l) Rc (3| + H.c, (38) 
3=1 

where the joints of the scattering center are L = 1, i? = 3 
and 

ff c = -« (|l) c (2| + |2) c (3| + |3) c (4| + |4) c (1| + H.c.) 
+z 7 i |2) c (2| - i72 |4) c (4| . (39) 

Note that here we consider a non-'PT'-symmetric model 
without losing the generality. The Bethe ansatz wave- 
function has the form 

—00 4 +00 

m = £ ^ \j) L + £ ^ \j) c + £ ^ \j) R , (40) 

J = -l 3 = 1 3 = 1 



where fj is in form of Eq. (|10|) . Taking k = 1, the explicit 
form of Schrodinger equations are 



— /_1 — /l2 — /14 = £71 1, 

-/l - h 2 - h 4 = £7i 3 , 

-hi -h 3 = (E- iji) h 2 , 
-hi - h 3 = (£ + ij 2 ) hi, 
E = — 2cosfc. 

the continuity of the wavefunctions demands 



(41) 



or equivalently in the explicit form as 

(H A ) mn = A (m| H VT \n) A = {H A )* nm , (33) 

(H B ) mn = B (m| ffpr = (H B y nm , (34) 
{H A B) mn =A (m\ H vr \n) B = zlm (V m ) <5 m _Ar lin ,(35) 
(H B A) mn =b (m\H VT \n) A = Urn (V n ) 8 m ,n-Ni . (36) 

Therefore, the cluster if-px ac ts as a Hermitian scattering 
center. This result is independent of the magnitudes of 
the hopping integrals and the potentials, also the reality 
of the spectrum of H-pt- 



hi = 1 + r, h 3 = t, 



(42) 



The corresponding transmission and reflection coeffi- 
cients have the form 



where 



c 



C cos k - 1 

i( sin fc ik 
t = —. -e , 

„-ik _ f ' 



1 1 



cos A; + £71 /2 cos /c — 172/2 
Straightforward algebra shows that 



(43) 
(44) 

(45) 



IV. ILLUSTRATION 

We consider a simple 4-site non-Hermitian scattering 
center which is illustrated schematically in Fig. HJa). 
The Hamiltonian of the whole system of Eq. ([T]) can be 



1 + ICI 2 - 2Re (C) cos fc + 21m (C) sin fc ' 

(46) 

which indicates the current is conserved when £ is real, 
i.e., 71 = 72. 
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Alternatively, taking the linear transformation 

|A) = (|2) C +|4) C )/V2, (47a) 
|B) = (|2) C -|4> C )/V2, (47b) 
the Hamiltonian He can be rewritten as 

H c = -V2k (\A) c (1| + \A) C (3| + H.c.) (48) 
+ . (7i-72) {A{ + ]B) m 

+ . (7i+72) (B| + |jB) 

as illustrated schematically in Fig. [2jb). It depicts a scat- 
tering configuration with single side coupling site which 
has been systematically studied in the Hermitian regime 
24]. Obviously, when 71 = 72 — 7, it is a simple ex- 
ample of a real Hermitian graph with parity-symmetry 
combining with the on-site PT-symmetric potentials and 
admits the current preserving. 

Accordingly, the transmission probability (coefficient) 
has the form 



T(k) 



sin 2 (2k) 



sin 2 (2fc) + (cos 2 (fc)- 7 2 /4)' 



(49) 



which has peculiar feature in contrast to that of Her- 
mitian scattering center. As a comparison, we write the 
transmission probability for the real side coupling by sub- 
stituting 7 with £7, i.e., 



r (k) 



sin 2 (2k) 



sin 2 (2k) + (cos 2 (k) + 7 2 /4) 



2 ■ 



(50) 



It can be observed that, (i) both of them have the com- 
mon total reflection points, T(tt/2) = T' (tt/2) = 0; (ii) 
T (k) = 1 at the resonance condition 7 = 2 |cosfc|, while 
T' (k) is always less than 1 within the whole range of 7. 



V. SUMMARY AND DISCUSSION 

We have proposed an anti- Hermitian coupled two Her- 
mitian graphs as the scattering center, which has been 
shown to act as a Hermitian graph, preserving the tradi- 
tional probability. This conclusion can be applied to the 
non-Hcrmitian scattering center which consists of pairs 
of PT-symmetric on-site potentials. This fact indicates 
the balanced gain and loss can result in the Hermitic- 
ity of the scattering center. Our results can give a good 
prediction for the transmission and reflection coefficients 
of linear waves scattered at the PT-symmetric defects 
in the experiment. The recent observation of breaking of 
PT symmetry in coupled optical waveguides [25|-[27j may 
pave the way to demonstrate the result presented in this 
paper. 

Finally, we would like to point that our conclusion can 
also apply to other type of non-Hermitian scattering cen- 
ter. For instance, we can select i? 7 , H a , H ia and H a fj 



being all Hermitian instead of real Hermitian in Eq. (|26l) , 
and we note that H-pj- is no longer 'PT-symmetric if the 
hoppings are not all real, with iJ 7 , H a , H la and H a p 
being Hermitian, we could also select H-pj- as 



H-, a 



H VT = H\ a H a + H & 



H 



10 

T 



H 



H a j3 

H„ — Ha 



(51) 



which is also in form of Eq. © after the transformation 
of Eq. (j2"5|) and exhibits the Hermitian behavior. 



VI. APPENDIX 

In this Appendix, we will prove the relation of Eq. 
2Tj) . For an incident plane wave with real energy E, we 



obtain from Eq. (|13|) and Ha 



H Al H B 



.ffjj that 



H A ~E -H AB 



H 



t 

AB 



Hb — E 



(52) 



Considering the block matrix A and A''' , when (Hb — E) 
is invertible, employing the Leibniz formula, we have 



det (A) = det (H B - E) det[(H A - E) 
-H AB (H B - E)- 1 (-H AB ) 



and also 



det (A f ) = det (H B - E) det[(H A - E) 



(53) 



(54) 



-(-H AB )(H B -E) X H 



AB\ 



Then we have 



det (A) = det (At) = [det (A T )] * = [det (A)]* , (55) 

i.e., det (A) is real. Such feature arises from the special 
structure of the matrix A in the form 



A C 
-Ct B 



(56) 



with A and B being Hermitian matrices. 

Matrix (i,j G [1, N A ]) is obtained from A by elim- 
inating its ith row and jth column, which has the form 



A' C 
V B 



and accordingly (Mji) has the form 



A' -C 
-V B 



(57) 



(58) 
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Here, A 1 is the matrix by eliminating the ith row and 
jth column from A, while C (T>') is the matrix by elimi- 
nating the ith row (jth column) from C (T>) . By similar 
procedure in obtaining Eqs. ([53")) and we have 

det(Afy) = det[(M ji ) t ], (59) 

and then 

det {Mij) = [det (Mji)}* . (60) 



Together with Eq. ([TSr). we have (A -1 ).. = (A" 1 )^, 
and yield Eq. (J22j). 
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